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Abstract. 

We consider singular solutions of a system of two cross-coupled Camassa-Holm (CCCH) 

equations. This CCCH system admits peakon solutions, but it is not in the two-component CH 
integrable hierarchy. The system is a pair of coupled Hamiltonian partial differential equations for 
two types of solutions on the real line, each of which separately possesses exp(— peakon solutions 
with a discontinuity in the first derivative at the peak. However, there are no self-interactions, so 
each of the two types of peakon solutions moves only under the induced velocity of the other type. 
We analyse the 'waltzing' solution behaviour of the cases with a single bound peakon pair (a peakon 
couple), as well as the over-taking collisions of peakon couples and the antisymmetric case of the 
head-on collision of a peakon couple and a peakon anti-couple. We then present numerical solutions 
of these collisions, which are inelastic because the waltzing peakon couples each possess an internal 
degree of freedom corresponding to their 'tempo' - that is, the period at which the two peakons of 
opposite type in the couple cycle around each other in phase space. Finally, we discuss compacton 
couple solutions of the cross-coupled Euler-Poincarc (CCEP) equations and illustrate the same 
types of collisions as for peakon couples, with triangular and parabolic compacton couples. We 
finish with a number of outstanding questions and challenges remaining for understanding couple 
dynamics of the CCCH and CCEP equations. 
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1. Introduction 

1.1. Purpose and scope of the paper 

In this paper we investigate solution properties of a new two-component Camassa-Holm (CH) 
system introduced in equations (13) and (14). As far as we know, this system is not completely 
integrable. However, the system does possess bound pairs of peakon solutions which exhibit 
interesting propagation dynamics involving both propagation and oscillation, while a single peakon 
must remain stationary (fixed in space). The oscillating and translating motion of the bound pairs 
of peakons as they propagate is reminiscent of the swirling motion of waltzing dancers, so we 
call these solutions peakon couples. The variational derivation of the equations, their geometry 
and the behaviour of the peakon couple propagation and collision interactions are studied in the 
paper. We also discuss compacton couple solutions for the more general case of cross-coupled 
Euler-Poincare (CCEP) equations in (13) and (14), and illustrate the same types of collisions as 
for peakon couples, with triangular and parabolic compacton couples. 

CH equation and generalizations The well-known CH equation may be written as [2, 3] 



where subscripts denote partial derivatives of the horizontal fluid velocity u{x,t) and the linear 
dispersion parameter a; is a real constant. The solution u{x, •) and its derivative Ux{x, •) are 
assumed to vanish sufficiently rapidly at spatial infinity, \im\x\^^{u,Ux){x, ■) —0. 

The integrable CH equation (1) describes the unidirectional propagation of shallow water 
waves over a flat bottom at one order in asymptotic expansion beyond the KdV equation 
[2, 3, 8, 9, 21, 22, 5]. The CH equation also governs axially symmetric waves in a hyperelastic rod 
[7]. In the case of no hnear dispersion, u — 0, the CH equation (1) possesses singular solutions in 
the form of peaked travelling waves called peakons 



Peakon properties and solution behaviour are reviewed in [11]. Recent developments about the 
CH equation may also be found in [18, 12] and the references therein. 

The scope of mathematical interpretations of CH may be gleaned by rewriting it in various 
equivalent forms, many of which have been rediscovered several times and each of which can be a 
point of departure for further investigation. For example, CH may be treated variously as: a fluid 
motion equation; a mathematical model of shallow water wave breaking; a transport equation 
for wave momentum by a fluid velocity related to it by inversion of the Helmholtz operator; a 
nonlocal characteristic equation; an Euler-Poincare equation describing geodesic motion on the 
diffeomorphism group with respect to the metric defined by the norm on the tangent space of 
vector fields; a Lie-Poisson Hamiltonian system describing coadjoint motion on the Bott-Virasoro 



Ut — Uxxt + 'iujUx + 3uUx — 




(1) 




(2) 
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Lie group; a bi-Hamiltonian system; a compatibility equation for a linear system of two equations 
in a Lax pair, etc. 

In particular, the CH equation on the real line follows from Hamilton's principle 6S = with an 
action S — J l{u)dt given by Lagrangian 

l[u) f u'^ + 2uu + ul dx (3) 



2. 

for solutions u{x, t) that vanish sufficiently rapidly at spatial infinity. The Euler-Poincare equation 
[16] 

51 

dtm = — ad*m = — {um)x — mUx with m — — u — u^x + i^, (4) 
then recovers the CH equation (1). 

The present paper emerged from the observation that complexifying the CH Lagrangian in (3) 
yields 

l{u, v) = - I {u + ivY + 2uj{u + iv) + {u^ + iv^Y dx 



~ 2 j ~'~ ^^■^ ~'~ ^'^^ ~'~ ~'~ ^^-^ dx j uv -\- UxVx + ojv dx (5) 
= lk{l)+i%l). 

Now the Euler-Poincare equation (4) for the real part of the Lagrangian 3?(/) simply yields two 
separate copies of the CH equation, one with ui and one without. This, of course, is of no further 
interest, because it gives nothing new. However, the Euler-Poincare equation for the imaginary 
part of the Lagrangian yields a new coupled system of two equations for u and v that we 
shall investigate in the remainder of the present paper. This system of cross-coupled CH equations 
(CCCH) is found to be (with details of the derivation given in Section 3) 

51 

dtm — — ad*m = — {vm)x — mVx with m :— — — u — Uxx + i^, (6) 

6v 

SI 

dtn — — ad*n = — {un)x — nux with n :— — — v — Vxx- (7) 

In these equations the momentum m (resp. n) is transported only by the opposite induced velocity 
V (resp. u). Hereafter, we ignore hnear dispersion by setting set a; = 0. This has the effect that: 
(i) the system is symmetric and (ii) it admits peakon solutions. Peakon solutions for the CCCH 
system (6) and (7) with ui — Q are singular solutions in which the values of the momenta m and 
n are supported with time-dependent weights on delta functions at points moving along the real 
line, carried by the appropriate velocity. 

The present paper focuses its attention on the dispersionless cross-coupled CH system (CCCH) 
comprising (6) and (7) with u — These equations govern two types of peakons. Each type of 
peakon is swept along by the induced flow of other peakons of the opposite type. We also consider 
solutions of the cross-coupled Euler-Poincare equations CCEP, 

dtm — — ad*m = — {vm)x — mVx with u — Kg * m, (8) 
dtu — — ad*n = — {un)x — nUx with v — Kg * n. (9) 
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For a convolution kernel Kg{-) with compact spatial support, these CCEP equations admit 
compacton solutionsi whose scattering interactions we shall compute. 

To finish establishing the notation, the next paragraph will quickly review the well-known 
behaviour of the CH equation by writing in several of its various forms and discussing its properties 
that will be relevant here, in parallel with the properties of the CCCH system (6) and (7). 
We should mention that there are other, integrable two and multi-component CH systems with 
application in fluid mechanics, see e.g. [23, 4, 19, 10, 6, 17, 20, 14, 13]. 



1.2. Parallel properties for CH and CCCH 

Momentum conservation in CH and CCCH. The CH equation (1) with a; = 

Ut — Uxxt + SwMa; — 2UxUxx — UUxxx = 0, (10) 

may be rewritten in the form of a momentum conservation law, as 

ut + dx(lu^ + p) =0, (11) 



^2 

with pressure P given by the convolution 

P — K * + 2^1^ with kernel K{x, y) = - exp(— |a; — y\)- (12) 

The kernel K for peakons is the Green's function for the ID Helmholtz operator, (1 — d^), and is 
also the shape of the peakon velocity profile in (2). 

When expressed in terms of the corresponding velocities u — K * m and v — K * n, with kernel 
K{x) given by the Green's function in equation (12), the cross-coupled CH (CCCH) equations (6) 
and (7) with ou — expand into 

Ut - Uxxt + VUx + 2VxU - 2VxUxx " VUxxx = 0, (13) 
Vt - Vxxt + UVx + 2UxV - 2UxVxx " UVxxx = 0, (14) 

As for CH, the total momentum of CCCH is conserved, namely 

dt {u-\-v)-\- dx{uv -\- K * {2uv -\- UxVx)) — 0. (15) 



Characteristic forms of CH and CCCH. The CH equation (10) may be interpreted as the 
condition for the 1-form density m dx"^ to be preserved (or 'frozen') in the flow of the characteristic 
velocity dx/dt — u{x{t),t), namely, 

^(jndx^^=0, along = u{x{t),t) = K * m. (16) 

Note that the relation between the velocity of the flow u and the property that it carries m is 
nonlocal. 

I The term 'compacton' was introduced in [24] to describe solutions of nonlinear evolutionary partial differential 
equations that propagate coherently, interact essentially elastically and have compact support. 



Two-component cross-coupled CH2 



5 



Likewise, for the CCCH equations, but with the cross velocities, we have 

= , along — = v{x{t),t) = K * n, (17) 

CLv 

)dx 
— 0, along = u{x{t),t) — K * m, (18) 

where again the characteristic velocities depend nonlocally on the quantities that they carry. These 
formulas imply that the signs of the momenta m and n are preserved along characteristics during 
the motion. In particular, if the initial conditions for either of the momenta m and n are everywhere 
positive, then the sign of that momentum is preserved throughout the subsequent evolution. 




Multi-peakon solutions of CH and CCCH. For CH, the velocity is related to the momentum by 
convolution with the kernel K in (12), which is the Greens function for the Helmholtz operator. 
When expressed in terms of the momentum, the multi-peakon solution of CH for a; = is a sum 
over delta functions, supported on a set of points moving on the real hue. That is, CH for a; = 
has solutions for velocity and momenta given, respectively, by 

N N 

uix, t)^Yl Pa(^)e-''~'"^*^' and m{x, i) = ^ Pait)5{x - g„(i)) . (19) 

a=l a=l 

As shown in [15], the CH peakon solution (19) is geometrically the cotangent-lift momentum map 
for the left action of the diffeomorphisms Diff(R) on a set of N points on the real line. The same 
momentum map applies for CCCH, resulting in the peakon solutions, with velocities 

M N 

^(^' ^) = 2 5^ """^^^ e-l^-'^''^*)! , v{x, i) = 2 E ^''(^) e-l""^^^*)' , (20) 

a=l 6=1 

and momenta supported on delta functions moving along the real line, 

M N 

m{x,t)^^ma{t)5{x-qa{t)), n{x,t) ^ ^nb{t) S{x - rb{t)) . (21) 

o=l 6=1 

The dynamics of the CCCH equations and their peakon solutions will be the focus of the remainder 
of the paper. 



1.3. Main content of the paper 

Section 2 studies the two-component solutions of the CCCH equations (13) and (14). This is done 
first in a numerical solution of their partial differential equations that establishes their characteristic 
wave behaviour in the classic dam-break problem for water waves. Then the study of their peakon 
solutions begins. 

Section 3 derives the system of cross-coupled CH (CCCH) equations in (6) and (7) from a 
variational principle on the direct- product space of C° vector fields X{M.) x X(M). 

Section 4 discusses A^-peakon solutions of the CCCH system. 
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Figure 1. Spacetime plots of the dam-break solution behaviour of the velocity variables, v (Panel 
a) and u (Panel b) of the cross-coupled equations (13, 14) for initial conditions (22) with tanh- 
squared profile and disturbance width 10 units. The subsequent flow is calculated numerically 
in a periodic domain of width 80 units. As in the classical Saint Venant shallow water problem, 
wave-like phenomena are created at the moment of release, though the waves propagate only in a 
rightward direction due to positive mean velocities. 

Section 5 studies the simplest possible case, M = N = 1. This is the peakon couple solution for 
CCCH. 

Section 6 discusses the periodic equilibrium solutions of the CCCH system. 

Section 7 considers the anti-symmetric collision of a CCCH peakon couple with its spatial reflection. 

Section 8 presents the results of numerical integrations of the interactions between the cross-coupled 
peakon couple and other coherent couple structures for CCCH. 

Section 9 studies compacton-couple solutions of the CCEP system (8) and (9). The dynamics of 
these compacton couples have some interesting features that depend on the shapes of their profiles. 
We study triangular and parabolic profiles. 

Section 11 recaps the main results of the paper and lists some of the outstanding problems 
remaining for future investigations of the CCCH and CCEP equations. 

2. Two-component solutions of the CCCH equations 

'Dam-break' equivalent problem for the cross-coupled equations. The dam-break problem involves 
a body of water of uniform depth, initially retained behind a barrier, in this case at x = ±5 in a 
periodic domain of width 80 units. When the barrier is suddenly removed at t = 0, the water flows 
downward and outward under gravity. The problem is to find the subsequent flow and determine 
the shape of the free surface. This question was addressed in the context of shallow-water theory, 
e.g., by Acheson [1]. The dam-break results for the cross-coupled system in (13)-(14) in Figure 1 
show evolution of the variables v (left panel) and u (right panel), arising from initial conditions 
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-20 20 -20 20 



(a) ^ (b) ^ 

Figure 2. Waterfall plots for the dam-break solution behaviour of the velocity variables v (Panel 

a) and (/ (Panel b) of the cross-coupled equations (13, 14) for the same data as Figure 1. Solutions 
at intervals of 3 time units are plotted sequentially up the page and the dots identify locations of 
the three highest local maxima at each time interval. 

Panel (a) shows the creation of a slow main excitation, a faster excitation and a few slower subsidiary 
excitations that branch off from it. 

Panel (b) shows a fast rightward moving excitation induced by the development of the first fast 
excitation in the left panel, and several slow subsidiary excitations that branch off from it. 

Both of these excitations are coherent structures that seem to retain their identities after the 
collision with each other in the middle of the figure. However, no compelling evidence about the 
integrability of the CCCH system can be drawn from this data. 

in a periodic domain, in which an initially localized disturbance with tanh-squared profile in the 
momentum variable n interacts with an initially constant mean flow in the independent cross- 
coupled velocity field, u. That is, 

n{x, 0) = [1 + tanh (x + 5) - tanh {x - 5)f , (22) 
m(x, 0) = u{x, 0) = 1. 

This dam-break initial condition first spawns a rapid leading pulse in velocity u and then produces 
a series of slower pulses in u, some of which are of larger amplitude than the first pulse. The 
rapid first pulse in u then passes again through the periodic domain, overtaking and colliding with 

the series of slower but larger pulses and undergoing a strong interaction indicated by a burst of 
amplitude in each collision. The u pulse and the v pulse track each other, but their variations have 
opposite phase. That is, a maximum in the u pulse corresponds to a minimum in the v pulse, and 
vice versa. This is an interesting scenario that we will see emerging when we study the coupled 
peakon solutions. 

3. Vciriational derivation of the CCCH and CCEP equations 

In this section we derive the CCCH and CCEP equations (6) and (7) from Hamilton's principle, 
by using the Euler-Poincare theory for symmetry reduction of right-invariant Lagrangians defined 
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on the tangent space of a Lie group. We begin by considering motion on the direct product of two 
copies of the smooth invertible maps with smooth inverses (diffeomeorphisms) acting on the real 
hne Diff(M) x Diff(R). We consider a reduced Lagrangian l{u, v) : X{R) x R that depends 

on a pair of smooth vector fields (u.v) G X(]R) = TgDiff (M)/Diff(M) that are right-invariant 
under Diff(]R). Applying the standard calculation of Hamilton's principle for this right-invariant 
Lagrangian yields, see e.g. [16, 18], 

l{u,v)dt^ f / ^,6u) + / ^,6v\ dt 

J a 



5u I \ 5v 

61 dr] , \ / 61 d^ , A , 

d 61 61 \^/d6l6l\ 

where one invokes rj = and C = for the variations in X(M) vanishing at the endpoints in time 
when integrating by parts and takes the continuum velocity vector fields {u, v) G X(]R) to vanish 
asymptotically in space. We have also used the formulas 

Su — ^ — adu77 , 5v — -J- — ad^^ , with rj, ^ E 3£(R) , (23) 

(JjL (Jjv 

for the variation of the right-invariant vectors field {u, v) and have taken the pairing 

(•,•): X(R) X X*(R) ^R 

as the pairing between elements of the Lie algebra X(R) and its dual 3£*(R), the space of smooth 
1-form densities. This pairing allows us to define the ad*-operation as 



where 



ad^^ := -[v, := -vC^, + ^v^ (25) 



is the adjoint Lie-algebra action X(]R) x X(]R) — >■ 3£(R), given by minus the Jacobi-Lie, or 
commutator bracket of smooth vector fields. In our case, (• , •) is the pairing between vector 
fields X(R) and 1-form densities X*(R), written as 

in which the integral is taken over the infinite real line. 

Substituting the reduced Legendre transformation u fi = 6i/6u and v ^ v = 61/ 6v produces the 
following pair of Euler-Poincare equations governing coadjoint motion for (/x, v) G X*(R) x X* 



dtjJ, -\- ad*// = , dtv -\- ad*i/ = for ji :— 6l/6u , u :— 6l/6v . (26) 
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We may evaluate the ad*-expressions from their definitions, so that 

(adV , ^) = , aduO = , + ^u^) = {{lJ'U)x + fJ-Ux, > 

and similarly for {n,u) — )■ (v^v). These expressions yield the Euler-Poincare equations on the real 
line for an arbitrary right-invariant Lagrangian. 

If the reduced Legendre transformation {u, v) {jj,, v) is invertible, then the Euler-Poincare 
equations (26) are equivalent to the (right) Lie-Poisson Hamiltonian equations: 



fi — -ad|^/5^// , 



(27) 



where the reduced Hamiltonian h and its variational derivatives are obtained from the reduced 
Legendre transformation by 

h{ijL, u) = (/X, u) + {v, v) — l{u, v). 

These equations are equivalent (via Lie-Poisson reduction and reconstruction) to Hamilton's 
equations on r*Difr x T*Difr relative to the Hamiltonian H : T*DiS x r*Difr R, obtained 
by right translating h from the identity element to other points via the right (diagonal) action of 
Diff on T'*Diff x T'*Diff . The Lie-Poisson equations may be written in the Poisson bracket form 



F{i,,u) = {F,h}, 



(28) 



where F : X* x 3£* — > R is an arbitrary smooth function and the bracket is the (right) Lie-Poisson 
bracket given by 



{F,G}{f,,u)^(fi 



'6F 


5G' 




'5F 5G' 


6fi ' 


6fi 


5v ' 5v 



(29) 



In the cases we shall consider, the Lagrangian l{u,v) will be bilinear, and its variations will 
cross-couple the momentum equations, as follows. An arbitrary norm on the tangent space 
of vector fields, = {u,Qu) , with Q a positive symmetric operator on u, induces cross- 
coupled equations through the application of the Euler-Poincare framework to the new Lagrangian 
Ig {u, v) — {u, Qv) — ( Qu, v) . In particular, the momenta in the Legendre transformation are given 
by, 

[I — = yv and v — —- — yu. 
6u 6v 

The corresponding Hamiltonian is given by 

hg{lJ^, u) ^ {Kg* ^ {Kg*u, ii) 

by symmetry of the kernel Kg. Consequently, the associated velocities are obtained from the 
convolutions 



6hg 

6n 



Kn * I' — u and 



Shg 
1^ 



Kg*fX 



Two-component cross-coupled CH2 



10 



in which the kernel Kg for the system is the Greens function for the invertible operator Q. The 
resulting Euler-Poincare or Lie-Poisson equations may both be expressed as 

dtGv + adlQv = , dtSu + adlQu = . (30) 

These are the equations of the CCEP system (8) and (9) for a general choice of the operator Q 
and its Greens function kernel Kg. They will be the CCCH equations (6) and (7) with u — 0, for 
the choice of the Helmholtz operator, Q — 1 — d^. 

Let us consider a few choices of the bihnear reduced Lagrangian i{u,v) whose Euler-Poincare 
equations (30) support peakon solutions on the real line: 

(i) When the reduced Lagrangian depends on only one vector field as £{u) = then the 
dispersionless CH equation, 

dtm + urrix + 2mUx = , with m = u — u^x, 

emerges as the dynamics of geodesic motion on the diffeomorphisms with respect to the 
norm IImUz/i = / u'^ + u'^dx. 

(ii) When the reduced Lagrangian is taken as i{u, v) = J uv+u^Vx dx, we find the following pair of 
coupled Euler-Poincare equations, in this case for the symmetric positive operator Q = l — d^, 
the Helmholtz operator, 

dtin + 2f^m -|- f m^. = , dtu + 2uxn + urix = , (31) 

with m = 6l/6v = u — u^x = Qu and n = 6l/6u = v — v^x = Qv. Because the momentum 
density (m) for one velocity vector field {v) is Lie-dragged by the vector field (u) for the other 
momentum density (n), we call this system the cross-coupled equations. When v — u, this 
coupled system restricts to the dispersionless CH equation. (Dispersion may be introduced 
by adding constant shifts to variables m — u — Uxx + ^^i and n — v — v^x + ^^2-) 

(iii) When {u,v) e C and v — the reduced Lagrangian becomes the complex norm 
£{u) ~ j {\u\'^ -\- \ux\'^) dx and we may interpret the coupled system of cross-coupled equations 
(31) as the complex version of the dispersionless CH equation. 

(iv) The system of cross-coupled equations (31) may be written in Lie-derivative form as 

{dt + JC^){m dx'^) = and {dt + £„) {n dx^) = , (32) 

with two characteristic velocities, u — K * m and v — K * n, in which the kernel K is the 
Green's function for the Helmholtz operator, defined in (12). The cross-coupled equations 
may also be written equivalently in a form similar to equation (18), as follows, 

— {mdx'^)—Q along — — v{x.,t) — K * n ., (33) 

^{ndx^) —0 along ^ — u{x,t) — K * m . (34) 

These equivalent forms of the cross-coupled equations imply they admit particle-like delta 
function solutions for both m and n, as discussed in the next section. 



The remainder of the paper is devoted to studying the solutions of the CCCH system (31) and the 
more general CCEP system (30). 
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We have aheady asserted that the solutions of the CCCH system (31) may be expressed in the 
form (21) for the two momenta and (20) for their corresponding velocities. The delta- function 
representation of these momentum solutions is also the cotangent lift momentum map for the 
direct product Diff(]R) x Diff(M) of diffeomorphisms acting on the real line IR by composition from 
the left [15]. 

By general principles for momentum maps and in parallel with the peakon dynamics for the CH 
equation, the 2M + 2N variables (^a, m^), a = 1, . . . , M , and (r;,, n^), h — 1, . . . ,N /va the peakon 
solutions for position (20) and momentum (21) of the CCCH equations (13, 14) are governed by 
Hamilton's canonical equations for the Hamiltonian function, 

M,N 

H=-Y, ma(^)r^6(t)e-l«"(*)-"''(*)l . (35) 

a,b=l 

These canonical equations comprise the evolution equations, 

BR 1 ^ 

Qait) = 1^ = ^ E n,{t)e-\^'^^'^-^^^'^\ = v{q,{t),t) , (36) 



drua 
dH 



=1 

M 



hit) = ^ = - E^«We-i««w-'-''Wi ^ u{n(t),t) , (37) 



dub 2 

a=l 

for the positions of the peakons, and 

N 



dH 1 ^ Qy 

ma{t) = - ^ = o"^a 5" nft sgn (g„ - rb)e-l^''(*)-^''(*)l = - m„ — , (38) 

OQa 2 f-' dx x=qa 

0=1 

dH 1 ^ du 
nb{t) = - ^ = --nbJ]m„sgn(g„-r6)e-l<'''W-'-''WI = -n,— _ , (39) 



x=rb 



a=l 

for their canonical momenta. In the case originally envisioned that the solutions would be complex, 
for V = iZ e C we have N — M and Ua — rha- 

Conserved quantities for these equations include the energy H and the total momentum Xla(^a + 
no). The solutions of these peakon equations will be studied in more detail in the sections that 
follow. 



5. The Peakon Dipole Solution 

In this section, we consider the peakon solutions of the CCCH system in the simplest possible 
case, M = N = 1. Dropping superfluous indices for this case, we introduce new canonical position 
variables X = {q + r) /2, Y = q — r , respectively mean position of the peaks and their separation 
distance. The evolution equations in terms of the new variables are 
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Thus we can define the behaviour of the exponential function of the absolute separation of the 
peaks, 

-el^l =sgn(y)^^-^, (40) 
where sgn {x) denotes the signum function 

1 for X > 



sgn [x) — 9{x) — 9 {—x) — 
Meanwhile 



1 for x < 0. 



m = -n = sgn (Y) —e-\^\ = sgn {Y) E, (41) 

where E=if|^^gis the (constant) value of the Hamiltonian, that is to say the total energy of the 
coupled pair. Differentiating (40) again with respect to time gives 

^ (gl^l) = -sgn^ {Y)E + S (Y) {n - mf . 
On integrating for a particular signature of yj^^Q = Fq 7^ 0, 



\Y\ nionoe ^^°k'^ (no - mo) 
e' ^ ^ sgn (Fo) ^ 

where 



mo = m|j^o , no = Yq = Y\^^^ . 

If mo and no have the same signature, then eventually we will have \Y\ =0, regardless of the value 
of |Fo|- If '^o and no are of different signature, then provided sgn (Fq) (^o ~ ^^o) < the same 
fate will occur. If the signatures differ and sgn (Iq) (no — mo) > then the particles immediately 
separate and continue to do so subject to the logarithmic rate of equation (42). In this case the 
amplitudes of both particles grow linearly according to equation (41). 

At the point of coUision when |y| = we require both H = \mn and M = m + n to be conserved. 
Thus the absolute value of m — n is also conserved. If the particles are assumed to 'tunnel' through 
each other then the momenta are conserved and the signature of the whole problem in Y switches, 
with the forcing on the amplitudes of the momenta reversed. Equivalently the particles can be 
assumed to collide elastically and to exchange both the amplitude and type of their momenta. In 
this interpretation, the signature of the forcings on the particles is set by the initial configuration. 
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When mo and uq share the same signature the half period of their 'waltzing' motion can be found 
by setting Yq — and looking for when e'^' attains unity, namely t — 2{no — nio)/nom,o. It will be 
noted that at this time 

fmo~no\ 
m\^^r, mo-no = mo + moTio I I = n-o, 

"»o"o \ moTT-o / 

and similarly 

SO that the two types of peakons do indeed exchange momentum amplitudes over a half cycle. The 
maximum separation of the particles occurs when the momenta are equal, 

no - mo no + mo 

, m^n^ , 

nomo 2 

at which point 

|F| = ln|l+(^°-^^°)'' 



4mono 

For two particles of arbitrary initial separation 1^, the full evolution of the 'particles' with respect 
to each other is described in terms of the conserved energy, 

E = ^monoe"''*^°', 
and total momentum 

M — mo + no 

in terms of the time variable, t, by equations 

mo + sgn (Fo) Et for t < Tq 

m{t)^{ mo + sgn (Yo) E (To - {t - for T2i<t< T^i+i 

no - sgn (Yo) E (To + {t - Ts^+i)) for T2,+i < i < T2, 

no - sgn (Yo) Et for t < To, 

n{t)^{ no - sgn (Yo) E (Tq + {t - T^i)) for T2, < i < T^i+i, 

mo + sgn {Yo) E (To - {t - Ts^+i)) for T2,+i < t < T2,, 



\l{t) 



In {-\Ee + |sgn {Yo) t (no - mo) + el^»l) for t < Tq, 

In (1 -\E{t- T,f + \{t- Ti) - 8E) for T, < i < T,+i, 
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sgn(r (t)) = 



sgn (lo) , for t < To, 

sgn(ro)(-l)\ for < t < Ti+i, 



where we have defined a modified signum function, 

1 for Yo > 0, 

sgn (lo) = { sgn {no - nio) for Yq = 0, 
-1 for Yo < 0. 

Here To is the time until the initial collision, where if > or sgn (Yo) {no — mo) < 0, this time 
is given by 



_ sgn (Yo) {no - mo) + VM-8E 
2E ' 

and To = oo otherwise. The subsequent Tj's in the peakon-peakon and antipeakon-antipeakon 
cases are the times of the subsequent collisions, 

VM^ - 8E 

Ti+i — Ti — - 



E 




(a) x-Xit) (b) 

Figure 3. Evolution of the velocity (Panel a) and momentum (Panel b) solutions of a coupled 
peakon pair. 

Panel (a) shows the profile of velocity fields of a peakon-peakon couple with initial parameters 
Too = 10, no — I, Yq — (solid lines). The dotted path indicates the evolution of the two peaks in 
the frame travelling at the particles' mean velocity, X = {q + r) /2. For these initial conditions, the 
total period of one orbit of the cycle is T = 3.6. Also shown are the velocity fields at subsequent 
time intervals t = 0.45 + 1.8fc for k E N (dashed lines). 

Panel (b) shows locus of the magnitude of momentum versus position for the two particles. The 
circles and crosses indicate the respective positions of each particle in phase space at time intervals 
of 1.44 units. The peakons are initially superposed at a; = 0, so that the leftmost circle and cross 
represent the same instant, and similarly for each couple continuing rightwards. 
Animations showing the time evolution of both these images are available as supplementary material 
with the online copy of this paper. 
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(a) time (b) time 

Figure 4. Plots of momentum (Panel a) and position (Panel b) versus time for the peakon-peakon 
couple with initial conditions tuq = 10, hq = 1, qq — tq — 0. 

Panel (a) illustrates that, as per equation (41) the growth and decay of momenta are continuous and 
piecewise linear with respect to time in the interval between collisions, while the total momentum, 
= m + n, is conserved through collisions. Also shown is the discontinuity (and symmetry) in the 
first derivative of the momentum at the points of collision, which here occur every 1.8 time units. 

Panel (b) illustrates the "waltzing" nature of the collision, with the identity of the leading particle 
alternating periodically, with the locus of the position of one particle obtainable from the other 
under a phase shift. Also plotted is the mean position X = {q + r) /2, showing the nonlinear 
relations in the propagation speed of the structure. 



The final step is to consider the mean fiow of the system. This is given by 

2ME + ^ ^ 



X 



X2_j2St-Sgn(yo){r).o-mo) 
2ME 



M'^-\lE{t-Ti)-^M'^-9.E\ 



I.e. 



^^^^-1 / 2i;t-sgn(^,)(no+mo) \ +sgn(y^)tanh-^ («) 



M 

X (T,) - Xo 



tanh 



-1 ( 2E{t-Ti)-VM^-8E 



t < Tr 



M 



+ tanh 



1 I VM^-SE 
M 







T,<t< T, 



Figs. 3-5 illustrate this analytic solution in the both the peakon-peakon {E > 0) and peakon- 
antipeakon {E < 0) cases. Of particular interest in the peakon-peakon case is the nature of 
the collisions that occur. These collisions propagate the waltzing peakon couple onwards as 
a coherent phenomenon, possessing (in the tunnelling representation) continuous values for the 
particle positions, velocities and momenta, but with a discontinuous first derivative in the rate of 
change of momentum. 

Note that over one entire cycle of the waltz, the mean speed of propagation is given by 



X 



^i + l ~ Ti Jrp 



Xdt 



mono 



2 \mo — no I 



tanh 



1 ( Imp - no\ 
mo + no 
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(a) x-Xit) (b) ^ 

Figure 5. Evolution of the velocity (Panel a) and momentum (Panel b) solutions of a peakon- 
antipeakon coupled pair. 

Panel (a) shows the velocity fields of anpeakon-antipeakon couple at time t = with mo = 10, 
riQ = —1, Yq = (solid lines). The dotted path indicates the locus of the peaks, in the frame 
travelling at the mean particle velocity, X {t) = (g + f)/2, as t ^ ±oo. Also shown are the 
velocities at i = 1. 

Panel (b) shows the loci of the (m, q) and (n, r) particles in a stationary frame of reference centred 
on the point of collision. As in Figure 3, the dots and crosses indicated the particles location in 
phase space at regular time intervals. Here the best point of comparison is from t = 0, when the 
particles are superposed at q ^ r ~ 0. 

Animations showing the time evolution of both these images are available as supplementary material 
with the online copy of this paper. 



Using THospital's rule we find that in the hmit of zero period oscillations 



limX= hm !!g/^ ^o + ^o-(^o-^o) 

To->0 mo^no 2 \ (mo + Uq) 



( 



1 



no 
2 ' 



This is exactly the velocity of a single peakon solution to the CH equations, to which the cross- 
coupled equations collapse under the same limit. 



Triangular compacton couples. Similarly, we may consider a compacton couple with a triangular 
kernel and corresponding Hamiltonian 

H = -mn(l — a\q — r\)9{l — a\q — r\) (43) 

where a is a constant. For initial conditions g(0) = r(0) = 0, n{0) = Uq, m(0) = mo, no > rriQ the 
solution is: 

rriQ + TT-oC"'^"*/^ ' no + moe^"*/^ 



m{t) = 
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Here <t <Ti — {4:/ac) In no/mo where c = mo + no is the conserved momentum. For the second 
half of the cycle the solution is 



cno cmo 



no + moe'=«(*-^i)/2 ' ^ ^ + ^^e-ca(t-Ti)/2 ' 

for Ti < t < 2Ti and the constant K = 2« In ^ + ^^ar^M represents the spatial shift during the 
first half cycle. The full cycle has a period 2ri = — In 

•' ac mo 

6. The Periodic Equilibrium Solutions 

While the 1+1 soliton has only one equilibrium solution, with m — n and I — 0, the periodic 
solution has two. The first being the CH equivalent state with particles super-imposed, and the 
second with m and n multipeakons which are a half-step out of phase with each other. In such 
cases, the Greens function solution is given in the interval —L<x<Lhy 



m{x, = ^ m„ {t) 5{x-qa (t)) , u{x, t) ^J^"^- (^) """"^ ^^4sinhL^° ^ 



a=l a=l 



n{x, t)^J2^, it) S{x-n it)) , v{x, t)^J2^, it) ^osh(L \x nit)\) ^ 



6=1 6=1 



and by the periodicity conditions, m{x + 2L) = m{x), n{x + 2L) = n{x), etc., otherwise. 
For the case M — N — 1, upon dropping the indices a and b and defining once again 

Y{x) — q{x) — r{x), 

we see (eventually) that 

^2 / L — \Y\\ 

— 2tan-^ f tanh ^ j = -2£;sgn^ (F) tanh (L - . 

From this relation, we obtain the quadrature 

d\Y\ /I \ Scosh(L-|y|) 

' '-' ln(cosh(L-|y|)) + C e SnRx^^ 



dt \2 

where as before 



2 
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and the constant of integration C is found from the initial condition 



d\Y\ 



dt 



,cosh(L-|lo| 
sgn [Yq) {no - mo) 



t=o ^ ^^^^ 



7. The 2+2 PecLkon-AntipecLkon solution 

We consider now the anti-symmetric collision of two couples of peakons in which, using the labelling 
pattern indicated in Figure 6 we have the following reduction of the system (36)-(39): 

mi = -777.2 = m, 771 = —n2 = 71, = —q2 = q, ri = -r2 = r. 

The reduced dynamical system in this case has the canonical Hamiltonian form. 



dH 



-m e 



-\<l-r\ _ f,-\q+r\ 



). 



dn 2 
dH 1 

777 = — -— = -mn (sgn(g — r)e~'^~''' — sgn(g + r)e~'^"'"''') , 

77 = — -— = 77777 (sgn(g — r)e^''^~^' + sgn(g + r)e~'^'''''') . 

ar 2 ^ V / / 

The Hamiltonian for this system is 

H = ^mn (e-l«-^l - e-l«+''l) 

and its canonical Poisson bracket arises from the symplectic form ou — dq A dm -\- dr Adn. 
We choose to consider the colliding case 

sgn (77) = sgn (777.) > 

and 

sgn {q) = sgn (r) < 0. 

Initially we suppose that g(0) = r(0) —\nk<0 where from the energy conservation 



2// 

A;= Wl- 



7?7o7lo 



7?7o = 777(0) and no — n{0). Then we assume that > q{t) > r{t) until after some time T 
again q{T) — r{T) and another 'exchange' of the positions of n and 777 takes place. Under these 
circumstances the system (44)-(44) for < i < Ti has the form 

q — —ne^ sinh q , f — —me^ sinh q , m — mne^ cosh q , n — mne^ sinh q 
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Figure 6. This is the head-on coUision set-up used in our analysis. It has a rightward traveUing 
peakon-peakon couple in the left half plane and an antisymmetric counterpart on the right. Each 
particle couple undergoes "waltzing" collisions before reaching, by symmetry, a singular collision at 
the point x = 0. 

An animation that shows the time evolution of such a peakon-antipeakon collision is available as 
supplementary material with the online copy of this paper. 

The system conserves the energy H = —mne^ sinh q and can be integrated, giving: 

Qit) = n ) i ) { , r(t) = n — , (44) 

(A; + 1)2 - (A; - l)2e^2noH 
m(t) = — , n(t) = no - Ht. 45 

The time to the next exchange, Ti can be determined from q{Ti) = r{Ti), i.e. from the equation 

, (k + l)-(k- l)e-"o'=^i , , , 

no - HT, = kno jy--^ yz pfT . 46 

(A; 1) + (A; - l)e-"o'=^i ^ ^ 

If Uq > mo this equation has a real root, which is less or equal to Uq/H. For < t < Ti n{t) is 
decreasing and m(t) is increasing. 

The equations for the time following the exchange, i.e. when > r{t) > q{t) are 
q = —ne'^ sinh r, r = —me'^ sinh r, 
771 = mne'^ sinh r, n = mne'' cosh q. 
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(a) 





(b) 

Figure 7. Evolution of the momentum of the peakon couple in the anti-symmetric headon collision 
in terms of position (Panel a) and time (Panel b). 

Panel (a) shows a log-log plot of momentum versus position for the peakon couple in the peakon- 
antipeakon collision, calculated from a numerical integration of the governing ODEs (44)-(44), 
(47)-(47). Here m (0) = l,n (0) = 2, g (0) = r (0) = -2. This plot illustrates the evolution over two 
complete cycles of the 'waltzing' peakon. while also showing the near exponential increase in the 
momentum amplitudes over each full cycle. 

Panel (b) shows momentum versus time for similar data, with m (0) = 1, n (0) = 2, q (0) — r (0) = 
—3. The last stage shows the rapid change from behaviour very similar to the pure peakon couple 
to near exponential growth as the particles in the peakon couple approach each other. In both 
cases the behaviour of the antipeakon couple is identical modulo a reversal of signature. 



with solutions 



nil- H{t-Ti) J ' V ('^i + 1) - ('^i - l)e-'"i'=i(*-^0 J ' 

m(t) =mi- H(t-Ti), n(t) = H ^^'^ ^ 7 ^^\7x ^^^^ .Tt! , 

where nii = m(Ti), ni = n{Ti) from (45) and (45) are the new (updated) initial conditions, and 
the updated value for the former constant k is now 



1. 



mini 



These solutions are valid until the time of the next collision, Ti < t < T2 when again q{T2) = r{T2) 
or, T2 is a solution a similar transcendental equation: 

mT T^ i. (fci + 1) - ih - i)e-^^^(^^-^^) 
mi - HiT, - Ti) = fei^i (^^^i)^(fe^_i),-^.MT.-T.) ' (47) 

The total time of the two exchanges, i.e. for a full cycle is T = T2. 

Thus the evolutionary process of the system that emerges from these solutions between the 
collisions is an iterative map illustrated in Figure 7. 
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-10 10 20 30 40 50 -10 10 20 30 40 50 



(a) ^ (b) ^ 

Figure 8. Spacetime plots showing an example of the evolution of the u (panel a) and v (panel 
b) velocity fields, under the interaction of a peakon couple with a free resting soliton. Here all 
particles are initially of momentum mi = TO2 ^ ni — A, with initial positions qi — 1, q2 — 20, 
ri = —1. Comparison of the panels shows the beating amplitude and phase locking of the waltzing 
couple, which is interrupted when the n particle interacts with the second m particle. 
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Figure 9. The space-time track of the particles presented in Figure 8. One waltzing partner deserts 
its old partner on interacting with another particle of opposite type, with a barely perceptible change 
in the period of oscillation (tempo) in the waltzing couple. 

Animations that show the time evolution of this image and the underlying peakons are available 
as supplementary material with the online copy of this paper. 

8. Other Collisions 

We now present the results of numerical integration of the interaction between a peakon couple 
and other coherent structures. The numerical code solves the ODEs for the particle momenta 
and peak locations using an implicit midpoint rule integrator and a collision detection routine to 
identify the points of discontinuity in the momentum forcing term. This method can be used to 
investigate the short term behaviour of the (M, A^) peakon equations for cases with M and/or 
greater than one, for which we do not yet have an analytic description. 



Two-component cross-coupled CH2 



22 



80 r 
70 - 
60 - 
50 - 
E40- 
30 - 
20 - 
10 - 

o[ 



(a) 



20 
X 



6.4 
5.6 
4,8 
4.0 
3.2 
2.4 
1.6 
0.8 
0.0 



(b) 




Figure 10. Spacetime plots showing an example of the evolution of the u (Panel a) and v (Panel 
b) velocity fields, under the interaction of two waltzing couples of peakons. Here mi = ni = 8, 
7712 = 772 = 4. gi = 92 = 14. Ti = 1, 7*2 = 16. The peakon couples interact and exchange some 
momentum leading two a speeding up of the rightmost couple. However the collision also affects 
the "internal oscillation" energy of the two couples of particles, changing the period of oscillation 
of both constructs. 




Figure 11. The space-time tracks of the waltzing couples of peakons from Figure 10. This makes 
it clearer that the two waltzing couples exchange momentum without ever actually overlapping, as 
well as illustrating the increased period of the lefthand couple and decreased period of the righthand 
one. Although not obvious, the post-exchange mean propagation speeds of the two couples also 
differ from those pre-coUision. 

Animations showing the time evolution of this image and the underlying peakons are available as 
supplementary material with the online copy of this paper. 



8.1. Interaction between a peakon couple and a solitary peakon 

This is the case with M = 2, N = 1, sgn (mi) = sgn(m2) = sgn(ni) in which a peakon couple 
interacts with a single peakon. The results of a typical integration are shown in Figs. 8-9, in 
terms of both the generated velocity fields and the tracks of the particles. Since the effect of one 
peakon on another falls exponentially with distance, a tightly bound peakon-peakon couple feels 
little influence from a distant solitary peakon, and by the same measure causes little change in the 
solitary peakon itself. However once the couple cycles sufficiently near at hand for the distance of 
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the solitary peakon from the couple to be comparable to the maximum separation distance of the 
couple, the particle feels sufficient effect to be swept into motion. The ordering condition on the 
line is sufficient to guarantee that the formerly solitary particle is never overtaken by the other 
member of its type. In all integrations performed, the solitary peakon is swept up into a new 
waltzing couple, with the former partner abandoned, however we have no definite proof that there 
does not exist a parameter regime in which the three particles may form a coherent propagating 
structure. 

8.2. Overtaking between two peakon couples 

This is the case with M = 2, = 2, sgn (m^) = sgn(nfe), in which two peakon couples interact. 
In general this requires the mean propagation speed of the trailing couple to be larger than that of 
the leading couple, so that the two structures interact in a manner equivalent to the CH overtaking 
collision. The results of a typical integration are shown in Figs. 10-11. As in the case of interaction 
with a solitary peakon, the particles behave in the manner of two pure cycling peakon-peakon 
couples, so long as the separation of the centres of the couples is significantly longer than the 
the maximum separation distance within either couple. Within the parameter space explored, the 
subsequent interaction is always at a distance, with no extra entanglement of the loops of the 
two tracks. The resulting exchange of momenta leads to an acceleration of the initially slower 
leading couple, with concomitant deceleration of the trailing couple. However, this exchange is 
not perfect, with additional energy partitioned into or out from the oscillatory motion, above that 
necessary to change the speed of progression. Prom the observed behaviour it can be conjectured 
that the long time solution of the {rii + n2,ni) peakon cross-coupled equations for ni,n2 G N is 
to form a chain of rii peakon-peakon couples, well ordered by their mean propagation speed, with 
n2 near non propagating solitary peakons remaining. Due to the complex nature of propagation 
speed formula it is not, however possible to guarantee well ordering in either the momenta or the 
separation distances involved. 

9. CCEP coupled compacton solutions 

An arbitrary norm on the tangent space of vector fields, \\u\\ — {u, Qu) , with Q a symmetric 
operator on u and its derivatives, induces cross-coupled equations through application of the 
CCCH framework to the new Lagrangian Ig {u, v) = {u, Qv) . In particular, the convolutions 
u = Kg * m and v = Kg * n, define the kernel Kg for the system, which inverts the momentum- 
velocity relation in the definitions of the momenta in the Legendre transformation. 



Slg 




m = — and n 
bv 



The numerical method, used in the previous section to integrate the Hamiltonian ODEs governing 
the peakon solutions of CCCH, was extended to solve the system of evolution equations for general 
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Figure 12. Plot shows the profile of the triangular and parabolic kernels, Ki and defined in 
equations (49) and (49), as well as the CCCH peakon profile, K^. The parabolic profile, K2, has 
discontinuities in derivatives at x = ±4 but is smooth across the origin. The triangular profile, Ki, 
has discontinuities in its derivative at all three points, x € [—4,0,4]. 

Hamiltonians of the form 

M,N 

^H=J2 m,{t)n,{t)Kg{\qa - rj). (48) 

a,fe=l 

We consider in particular the two compactly supported kernels, 




that are, respectively, triangular and parabolic. The profiles of these two compacton kernels are 
displayed in Figure 12, in comparison with the peakon profile, ^exp(— |a;|). 

Much of the analysis of Section 5 generalises immediately to singular solutions of cross-coupled 
Euler-Poincare (CCPE) equations with compacton kernels. When the centres of two solitary m 
and n particles lie within each other's domains of support, the particles propagate together with 
the same form of cyclical motion for CCPE compactons as exhibited by the CCCH peakons. For 
example, numerical integrations shows the same results as for the peakon in the interaction of a 
cyclically propagating compacton couple with a solitary, stationary compacton, for each of the 
kernels, Ki and These cases are essentially the same as in Figure 9 for identical initial particle 
weights and separations. 



A new phenomenon for compacton collisions. The collision problem for parabolic compaction 
couples with the K2 profile does show an additional new property of the compacton systems. 
Namely, pair interactions may force the oscillations past the half- width of the compacton's support, 
and thereby separate a waltzing compacton couple into two solitary compactons. The role of this 
phenomenon in multi-compacton interactions is explored in order of increasing complexity in Figs. 
(13), (14) and (17). 



Two-component cross-coupled CH2 



25 




Figure 13. Panel (a) depicts spacetime plots for the pairwise overtaking problem equivalent to 
Figure 11, but with the parabolic compacton profiles in equation (49). The interaction is sufficiently 
violent that the trailing n (blue) particle gains enough extra velocity to escape the region of support 
of its partner (red), breaking the cyclic red-blue couple into two stationary particles. 
Panel (b) shows the positions and weights of the particles shortly after the collision, as well as 
the resultant velocity profiles. This clearly shows that the separation of the particles exceeds the 
compacton width. 

This behaviour was also found for other compacton kernels (such as triangular compactons) in 
regimes with sufficiently energetic compacton couples. Overtaking collisions of cross-coupled 
compactons can thus influence the form of the long time solutions. This behaviour is not seen 
for systems of compacton equations that contain only a self interaction term. 

Animations that show the time evolution of this image are available as supplementary material 
with the online copy of this paper. 



10. Crossflow equations with more than two species 



As a final example of the application of of this framework, we consider equations for multiple (i.e. 
more than two) species that cannot sustain motion without being in each other's presence. For a 
vector of velocity vector fields, u = (^u^^\u^'^\ . . . define a Lagrangian 

liu) = ^iu,iA) Qu), 

with A a real symmetric matrix of determinant ±1, and Q an arbitrary differential operator, 
applied elementwise. The resulting Euler-Poincare equations are 

for a vector of momenta, m := ^ = AQu. The inversion of this Legendre transformation is 

u = A-\K *m) = K * A-^m, (51) 
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Figure 14. Spacetime paths for the problem of sequential pairwise overtaking collisions are 
plotted for three initially rightward- moving parabolic compacton couples, kernel K2, arranged 
from left to right in order of decreasing propagation speed and one stationary parabolic compacton 
monopole. Initial conditions are given in form (species, weight, position) by (m,8, 0), (n,0.8, 0), 
(m, 1.6, 8), (n, 1.6, 9.8), (m, 0.8, 24), (n, 0.8, 24.8) and (m, 1.6, 48). During the overtaking collisions, 
the left-right positional order is preserved separately for each species. 
The first four collisions shown can be classified into three classes: 

(i) In the first collision, near t = 8, x = 10 the two colliding couples remain intact, but suffer 
sudden changes in their direction of propagation and in the periods of the cycles in both 
couples. 

(ii) The second collision at t = 22, x = 30, involves the destruction of the faster (dashed) couple, 
resulting in two unbound monopoles, one of each species, the second of which (dashed blue) 
is ejected only after propagating for a few cycles as a compound state involving two blues (m) 
and one red (n). The collisions at later times all involve exchanges of partners. 

(iii) The penultimate collision again involves formation of a compound state of two blues and one 
red. The compound state lasts for only two cycles in this case before the exchange is complete 
and the leftmost one (solid blue) is ejected. 

(iv) The last collision involves the stationary monopole that was ejected in an earlier exchange 
collision. It exchanges again and propagates away in a couple. Thus, the initially rightmost 
positive compacton couples propagate away rightward. 

Animations that show the time evolution of this image are available as supplementary material 
with the online copy of this paper. 



where the convolution is also to be applied elementwise. The induced Hamiltonian for this system 
is thus 

H =- (m, K * A~^m) , 
2 
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with Hamiltonian motion equations 

5H 

dtvn + idxTn + mdx) — = 0. (52) 

om 

An absence of self-interaction terms in the Lagrangian representation (50) requires that the matrix 
A have aU entries on its leading diagonal vanish. This is the necessary and sufficient sufficient to 
ensure that no non-linear terms in a single species appear in the governing PDEs for the velocities, 
meaning initial conditions with non-vanishing velocity of a single species are steady. 

Meanwhile, an absence of self interaction terms in the Hamiltonian representation (52) requires 
similarly that the inverse matrix have all the leading diagonal vanish, so that the the evolution 
of each element of m is instantaneously independent of its own value. 

For n — 2 the only matrix satisfying our conditions on symmetry, determinant and leading diagonal 
is 




which is unitary and generates generalizations of the CCCH. Thus, for this case of n = 2 the 
CCCH are the canonical equations without self-interaction terms. For n > 2 there may exist 
systems which are non-self interacting in only one of the two representations. 

As a concrete example, consider the case for n = 3 and Euler-Poincare Lagrangian 

3 



i.e. the coupling matrix for the Euler-Poincare case 

[A^i^ = ^ (1 - S,j) (53) 
has only cross-interaction terms. The choice of coupling matrix [A^^] gives momenta 



(1) SIep a d 
m> > :— ^ — 1 



m(2)+^,(3)' 



(54) 



and cyclic permutations of 1,2,3. The velocities then relate to the momenta by 

= (55) 



and cyclic permutations of the labels 1,2,3. The EP multi-species cross-coupled system has 
Hamiltonian 

Hep = ^ f ^^^^^ * ^^''^ ~ Yl ^^^^^ * ^^^^dx, (56) 
j^k=i j=i 
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in which the last term involves self-interaction. 

Another, complementary, cross-coupled system exists, in which the Hamiltonian representation has 
no self-interaction terms. This complementary cross-coupled system is a Lie-Poisson Hamiltonian 
system generated by 

i^LP = ^ / m^^^K * m^^Ux (57) 
with velocities given by, cf. equation (55), 

"U7,(2) _^ ^(3) 



(58) 



and cyclic permutations of the labels 1,2,3. This formula implies the following Lagrangian for the 
LP multi-species cross-coupled system. 



3 



j^k=l j=l 
and thus a different coupling matrix from (53), namely, 

[^'n^.= [^Ep],, = (l-25..) (59) 

Figures 15-16 show the results of integrating the singular solutions of these two different equation 
sets for the same initial conditions (in momentum space). The two equation sets have markedly 
different behaviour because of their differences in interactions. In particular, the self-interaction 
term in the EP Hamiltonian (56) induces a leftward drift for positively weighted particles and 
thereby allows for complexly coupled three-way interactions. In contrast, the three-species problem 
for the LP cross-coupled equations has fewer differences from CCCH for two species. 

As more species are added, the scope for more complex interesting interactions grows. For example. 
Figure 17 shows the result of a run with parabolic compactons of twelve separate species in the 
relevant Hamiltonian cross-flow system. This shows both the complex entanglement of the various 
particles and the trivial classicifaction of many interaction features into those already described in 
the two species compacton case. 

The primary stable features for the multi-species collisions in Figure 17 are the same as for two 
species collisions in Figure 14, in that both include solitary particles and waltzing pairs. However, 
the transient collision behaviour in the multi-species case can be significantly more complex, 
possessing, for example, fully entangled three-way collisions, particle tunnelling and "hop-scotch", 
in which a solitary particle is bled of its momentum by a coupled pair, then abandoned. 
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Figure 15. Initial conditions for three species cross-coupled equations in the Lagrangian (Panel 
a) and Hamiltonian representations (Panel b). In each case the locations and weights of the 
momentum particles, m^*', are identical. However, the resultant velocity fields, u^^\ which carry 
each of the particles differ markedly. Comparision of the panels shows the impact of the strong 
negative self-interaction term in the viscinity of each particle in the Lagrangian representation. 




(a) ^ (b) ^ 

Figure 16. Spacetime plots for three species cross-coupled equations in the Lagrangian (Panel a) 
and Hamiltonian representations (Panel b) with the initial conditions shown in Figure 15. 
Panel (a) shows the great significance of the negative momentum self-interaction term which appears 
in equation (54). All particles experience a leftward drift, proportional to their own magnitude. 
Moreover the fully coupled nature of the momentum equations allows. 

The Hamiltonian representation in Panel (b) behaves much more in the manner of the CCCH 
solutions. In fact for these initial conditions the solution is very similar to one in which the to*-^' 
and m^'^^ particles are of the same species. 
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X 

Figure 17. A space-time plot is shown of the evolution starting from initial conditions consisting 
of 12 separate species in the LP cross-coupled representation, whose Hamiltonian, as in (57), has no 
self-interaction terms. The initial conditions are chosen to form five parabolic compacton couples 
and two single parabolic compactons. All of the compactons are positive, so all of the couples 
propagate rightward. The same three basic forms of collision appear as in the previous figures for 
the CCCH equations however the particles no longer discriminate in their exchange of partners. 
The final space-time plot is reminiscent of the map of the London Underground! 
Animations that show the time evolution of this image are available as supplementary material 
with the online copy of this paper. 

11. Conclusions/Discussion 

Explicit outstanding problems. Several outstanding problems were identified explicitly in the 
course of this work. These include the following tasks. 

(i) Study and classify the full set of interactions of the peakons in (20) for the CCCH system 
(31). Determine whether the generalizations proposed in systems (31) are integrable, or even 
whether the (weak) solutions exist and are unique. 

Even if the CCCH equations are not completely integrable, the generalisation they represent 
is interesting for its singular solutions and geometrical interpretation. In this regard, we 
may remark about additional symmetry and complete integrability of the corresponding 
coupled rigid body dynamics. Namely, if the problem specified here for a right-invariant 
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Lagrangian on X(R) x 3£(R) had been expressed instead on 50 (3) x 50 (3) for a left-invariant 
Lagrangian, the result would have been interpretable as the dynamics of a cross-coupled 
system of two rigid bodies. This dynamics would have been expressible in Lie-Poisson form 
and it would have been completely integrable as a Hamiltonian system, only provided an 
additional 5*0(2) symmetry were present, as for the case of axisymmetric coupled rigid bodies. 
In the present case, there seems to be no extra symmetry that would make the corresponding 
CCCH Hamiltonian system (31) completely integrable. 

(ii) Explore the dynamical systems properties of the (M, A^)-peakon solutions (20) of the CCCH 
equations (31). For example, determine the results of colliding a compacton couple with a 
single compacton of opposite sign. 

(iii) Characterise the creation of peakons and peakon couples as a process of blow-up in finite 
time. We have sought an analog of the steepening lemma for the CH equation in [2] to 
explain how peakon couples arise from smooth initial conditions for the CCCH equations, 
but so far this result has eluded us. 

(iv) Discuss the solution behaviour of other cross-coupled equations, such as cross-coupled Burgers 
(CCB), cross-coupled Korteweg de Vries (CCKdV) and cross-coupled 6-equations. This is a 
wide open problem. 
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